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Summary. — In the seconds after core collapse and explosion, a thermal neutrino- 
driven wind emerges from the cooling, deleptonizing newly-born neutron star. If 
the neutron star has a large-scale magnetar-strength surface magnetic field and 
millisecond rotation period, then the wind is driven primarily by magneto-centrifugal 
slinging, and only secondarily by neutrino interactions. The strong magnetic field 
forces the wind to corotate with the stellar surface and the neutron star's rotational 
energy is efficiently extracted. As the neutron star cools, and the wind becomes 
increasingly magnetically-dominated, the outflow becomes relativistic. Here I review 
the millisecond magnetar model for long-duration gamma ray bursts and explore 
some of the basic physics of neutrino-magnetocentrifugal winds. I further speculate 
on some issues of collimation and geometry in the millisecond magnetar model. 



1. Introduction 

A successful core-collapse supernova (SN) leaves behind a hot deleptonizing proto- 
neutron star (PNS) that cools and contracts on its Kelvin-Helmholtz cooling timescale 
(tkh ~ 10 s), radiating its gravitational binding energy (~ 10 53 ergs) in neutrinos pp. 
A small fraction of these neutrinos deposit their energy in the tenuous atmosphere of 
the PNS through the interactions v e n — > pe~ , D e p — > ne + , vv — > e + e~ , and v(eT e + ) — > 
u'(e~'e +/ ). Inverse processes provide cooling. Net neutrino heating drives a thermal wind 
that emerges into the post-supernova shock environment pj]- For typical non-rotating 
non-magnetic (NRNM) neutron stars, the total kinetic energy of the wind over tkh is of 
order ~ 10 48 ergs and, hence, the addition to the asymptotic SN energetics is small on 
the scale of the canonical SN energy, E SN ~ 10 51 ergs 013 E]. 

Magnetars — a class of neutron stars with surface magnetic field strengths of B ~ 
10 15 G — are thought to be born with millisecond rotation periods, their intense fields 
having been generated by an efficient dynamo |SI • Millisecond rotation periods 
imply a reservoir of rotational energy that is large on the scale of -Esn: -E-Rot ~ 2 x 
10 52 ergs M 1A Rl w P^ 2 , where M 1A = M/1.4 M , R Vw = R u /W km, and Pi is the 
spin period in units of 1 millisecond (ms). Stellar progenitors that produce millisecond 
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magnetars (MSMs) must have iron cores rotating with periods of ^ 10 s at the moment of 
collapse [5]. The character of collapse, bounce, and explosion can be modified by rotation. 
For example, rotational support leads to characteristically lower neutrino luminosity 
(L u ) and longer tkh jlOl 1 1 1 j . In addition, a fraction of the gravitational binding of 
rotating collapse may be stored in shear energy and tapped by viscous processes; for 
parameters appropriate to MSM birth this effect can yield explosions in models that 
would otherwise fail Although rotation may be important during collapse and 

explosion, and although small-scale magnetic fields may be important in providing a 
viscosity capable of tapping the shear energy generated during collapse, even the magnetic 
energy density associated with a dipole field strength of 10 15 G is small with respect to 
the thermal pressure exterior to any PNS in the first ^1 — 2 seconds after collapse. 
Therefore, similar to NRNM PNSs, we expect that the wind that accompanies MSM 
cooling is driven by neutrino heating at early times. However, as L v decreases and the 
thermal pressure exterior to the MSM decreases, the region exterior to the MSM must 
become magnetically-dominated. 

The strong magnetic field forces the matter composing the outflow into near corotation 
with the stellar surface out to ~ Ra, the Alfven point, where the magnetic energy density 
equals the kinetic energy density of the outflow. For P ~ 1 ms, if Ra > 15 km, then 
the wind is driven primarily by magneto-centrifugal slinging; neutrino heating becomes 
relatively unimportant in determining the asymptotic wind velocity. Rotational energy is 
transferred directly from the MSM to the wind and this provides an efficient mechanism 
for spindown ^] • The spin period of the MSM e-folds on the spindown timescale 
tj w (2/5)(M/M)(R„/Ra) 2 , where M is the mass loss rate. Because -ERot ^> Esn, just 
one e-folding of is sufficient to modify the dynamics of the SN remnant significantly. If 
tj is small with respect to the time for the SN shock to traverse the progenitor (~ tens 
of seconds for type-lb, -Ic progenitors) we also expect this extra energy source to modify 
the SN nucleosynthesis As the MSM cools and the outflow becomes increasingly 

magnetically-dominated, Ra increases. It cannot do so indefinitely. Ra approaches the 
radius of the light cylinder i?L = c/il ~ 48Pi km asymptotically. As it does so, the flow 
becomes increasingly relativistic. This is the transition between non-relativistic mass- 
loaded outflow and relativistic Poynting-flux dominated neutron star wind. All neutron 
stars, regardless of their initial spin period and magnetic field strength, go through this 
transition. MSMs are interesting because this transition occurs at high wind kinetic 
luminosity. Because E^ ot is large and the spindown timescale is short, and because 
the velocity of the wind must eventually become relativistic, these objects are a natural 
candidate for the central engines of long-duration gamma ray bursts (GRBs) |15U16llT2*] . 

2. Proto-Magnetar Spindown 

In this section I summarize the results of Ref. [12] • See that work for more details. 

Angular momentum conservation implies that J = d/dt(ICl) = —MC, where C is 
the specific angular momentum carried by the wind and / is the moment of inertia. 
In the classic model for solar spindown constructed by Ref. t ne wind problem is 
treated in one spatial dimension and in the equatorial plane. Consideration of the az- 
imuthal momentum equation together with Faraday's law gives C = R^fl. To estimate 
the angular momentum loss rate and, thereby, the wind luminosity and the spindown 
timescale, we must first estimate Ra- The location of the Alfven point depends on the 
radial dependence of the poloidal magnetic field. Because models with purely monopolc 
fields generally over-estimate spindown and models with pure dipole fields under-estimate 
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spindown, we parameterize B r — Bo(R u /r) n , where 2 < r\ <J 3. Taking B 2 /8tt ~ pv 2 /2 
at i?A, P = M /Airr 2 v r , and assuming that va = v t (Ra) ~ v^(Ra) ~ i?A^, we find that 

(1) Jfc = B^ 2 "" 1 ) ^/(2„-i) (mjj)-!/^,-^ 

where wa is the radial Alfvcn speed, tv is the radial velocity, is the azimuthal velocity, 
p is the mass density, and B r is the radial magnetic field. Equation Q assumes that 
Ra& S> Uy, where ti„ is the asymptotic wind velocity in a NRNM outflow (y v ^ 3 x 10 9 
cm s _1 ). The absolute value of the rotational energy loss rate can be written as 

(2) E NR = j B4/(2')-l) i? 4r ; /(2r ) -l) j^{2r,-3)/{2r,-l) Q(4t7-4)/(2tj-l) _ 



For parameters appropriate to MSMs 

(3) ££r * 1.5 x 10 51 < 3 < 3 My 3 3 Pr 4/3 ergs s" 1 

(4) £^= 3 - 4.5 x ltf°Bgl%/*M%*P{*'* ergs s" 1 . 



The subscript 'NR' is added to emphasize that when the flow is non-relativistic, E 
depends explicitly on M. 

The non-relativistic scalings for the energy loss rate can be compared with those in 
the relativistic regime. As Ra becomes close to i?L as M decreases during the cooling 
epoch, va approaches c and the flow becomes relativistic. At _Rl, the ratio of magnetic 
to kinetic energy density is 



(5) 



B 2 



4irpc 2 



Rl 



and the energy loss rate is 



(6) 



E R = -TMc 2 = -BlR 2 jm 2l i~ 2 e'- 2n . 



For r\ = 3 the classical "vacuum dipole" limit is obtained; E R (rj = 3) = i?Qi?^f2 4 c _3 .( 1 ) 
For the monopole case with rj = 2, E R is larger than the dipole limit by a factor of 
c 2 /(f2i?„) 2 — a factor of ~ 23 for a 10 km MSM with a 1 ms spin period. 

The non-relativistic spindown rate is larger than the relativistic spindown rate as a 
result of mass-loading. To see this explicitly, note that 



(7) 
(8) 



E^ 2 /E R = 



r;=3 
R 



IB^R-^M^P?* 

U15 via —o 1 

8B- 6/5 R-^M% 5 Pl 2/5 . 

U15 via —A i 



For rj = 2 we see that the ratio is approximately unity, reflecting the fact that for 
the parameters chosen Ra ~ Rl- For slower spin periods, the ratio increases, but 
not dramatically. In contrast, the ratio of -Enr to E R for the dipole case (rj = 3) 



( x ) In the true "vacuum dipole" limit Eb, has a term sin a/6, where a is the angle between the 
spin axis and the magnetic dipole axis. 
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is large for MSMs and it has a strong dependence on P. Thus, for a magnetar born 
with a 10 ms spin period a naive application of the "vacuum dipole" formula underes- 
timates the magnitude of the rotational energy loss rate by a factor of ~ 2000. More 
detailed calculations reveal that when P ~ 10 ms, for modest neutrino luminosities, 
M is probably closer to ~ 10~ 5 M Q s _1 so that the ratio of £nr to Er, is closer to 
~ 120 than 2000. Even so, -Enr / E^ is very large and an application of the relativistic 
formula in an epoch when M is large is incorrect. Of course, if the field structure is 
not purely dipolar (that is, rj < 3), then the discrepancy between the "vacuum dipole" 
spindown approximation and the "true" spindown rate becomes even larger. In this case 
we would compare E^/E^ 3 ~ 25B~ i l /3 R~ 1 1 0/3 M 1 _ /3 P* /3 , a factor of 25 for a MSM. 
For a magnetar with a 10 ms spin period and lower mass loss rate the ratio becomes 
^nrV^r -3 — 2500B^/ 3 RZ^^ 3 M 1 /^ Pf^ 3 . These arguments serve to underscore the 
fact that in the very early stages of proto-magnetar cooling, the multi-dimensional struc- 
ture of the wind/magnetic field interaction must be solved consistently (to determine the 
effective rj) and that inferences about the "initial" spin period of magnetars should not 
be based on an application of the vacuum dipole approximation when M is large. 

For fixed Bq, i?„, and Q the transition from non-relativistic (eq. |3]) to relativistic 
(eq. jHj) outflow occurs when Ra ~ Rl and this point in time corresponds to a critical 
mass loss rate M C rit = B 2 R 2ll Vl 2,1 ^ 2 c 1 ^ 2ri , which scales with P~ 2 and P~ 4 for rj = 2 and 
77 = 3, respectively. Because M scales with the neutrino luminosity and because — to 
first approximation the luminosity is a monotonically decreasing function of time, this 
scaling of M C rit with P implies that the wind is non-relativistic for a larger fraction of 
the cooling time for longer rotation periods: for 77 = 3 and P = 10 ms M C rit — 3 x 10~ 9 
M Q s _1 . Such a low M may correspond to a time several tens of seconds after collapse. 



3. Are Millisecond Proto-Magnetars GRB Central Engines? 

The spindown timescale Q/fi in the non-relativistic limit is 

(9) r- 2 - 30 s M, 4 MZl 1 * R^J 3 B^* P^ 3 

(10) T^l - 96 s M 1A MZl ,b R^J 5 B^ 4 / 5 P^ 



For a MSM with 2 x 10 52 ergs of rotational energy we need only wait a fraction of tj to 
extract an amount of energy comparable to the supernova energy, ~ 10 51 ergs. Because 
tj (or ~ Tj/5) is comparable to tkh we infer that significant spindown may occur during 
the cooling epoch. As implied by the discussion of -Enr above, for larger initial spin 
periods the spindown timescales decrease. In the relativistic limit (Ra ~ -Rl) 

(11) t-jT 2 - 34 s M 1A R~l B£ 

(12) r^ 3 ~ 760 s Mx.4 R£ B£ P? 



Although Tj~ 2 does not depend on Q explicitly, it may have an implicit dependence on 
f2 if the magnetic field is generated by a dynamo [HIE]. 

With these timescales in hand we can (in a rudimentary way) attempt to assess the 
MSM GRB mechanism. A more detailed assessment must await a consistent multi-d 
MHD solution. From eqs. ©, ©, ©, and <|10Jl we see that in the non-relativistic limit, 
regardless of 77, the amount of energy extractable on ~ 10 — 100 second timescales is 
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in the range appropriate to GRBs. For rj < 3 these conclusions are stronger. Because 
> 10 51 ergs can be extracted on a timescale shorter than or comparable to the timescale 
for the SN shock to traverse the progenitor, we expect that the wind may significantly 
affect the dynamics of the remnant and the 56 Ni yield. In this way it may be possible 
to generate hyper-energetic or 1998bw-like SNe during MSM birth |12U18j . The inferred 
energetics and 56 Ni yield of SN2003dh and SN1998bw put strong constraints on any GRB 
mechanism. In the collapsar model a disk wind is thought to generate the 56 Ni required 
to power the SN lightcurve ^3 |20j ■ In the millisecond magnetar model, the energetic 
wind shocks the material already processed by the supernova shock, perhaps generating 
the large inferred 56 Ni yields; we are currently investigating the timing of this scenario. 

Although it seems possible that MSM winds may generate energetic winds at early 
times, the flow during this mass-loaded wind phase — at least, on average — is not 
relativistic. It is possible that a strong latitudinal dependence to the mass loss rate 
may yield relativistic asymptotic velocities for matter emerging from mid-latitudes even 
when our estimates would indicate i?A < -Rl, but such a speculation must be tested 
against realistic multi-d models. It is also possible that large temporal variations in M 
could cause the wind to alternate rapidly between non-relativistic and relativistic. Strong 
variations in the mass loading could be caused by shearing of large-scale closed magnetic 
loops on the surface of the fully convective MSM core |16j. 

As the flow becomes increasingly relativistic (on average), we see from equation Ijl2|l 
that if the relativistic dipole limit strictly obtains, then the spindown timescale is long 
for B ~ 10 15 G and P ~ 1 ms. Although relativistic spindown with rj = 3 will affect the 
asymptotic remnant dynamics by injecting energy over a long timescale, it will probably 
not generate a GRB with duration ~ 30 seconds and energy ~ 10 51 ergs. Based on 
the scalings derived here, higher magnetic field strength, shorter spin period, or 77 < 
3 is probably required in order for MSM spindown to power GRBs. If one of these 
possibilities obtains, then from the estimates in Ref. ^] we find that essentially all of 
the magnetic energy at Rl must be transferred to the wind in order for the flow to obtain 
high asymptotic Lorcntz factor with large enough E. This is presumably accomplished 
by magnetic dissipation |21l I22j . 

4. Emergence, Geometry, & Collimation 

If a relativistic outflow with the requisite energy to power a GRB can be generated 
by a MSM, it must emerge from the massive stellar progenitor. The highly energetic 
non-relativistic wind, which precedes the relativistic outflow, will likely be collimated by 
hoop stress and will therefore shape the cavity into which the relativistic wind emerges. 
Because the non-relativistic wind carries little mass in comparison with the overlying star, 
the relativistic outflow is not additionally hindered in its escape from the progenitor by 
the preceding slow wind. Hence, if the relativistic outflow can be collimated, then the 
dynamics of its emergence from the progenitor should be qualitatively similar to models 
of collapsar jets escaping Type-Ibc progenitors |23|. 

One important possible objection to the MSM mechanism for GRBs is that it is 
difficult to collimate relativistic Poynting-flux dominated outflows Observational 
evidence for collimation in GRBs is abundant and so, at face-value, this would seem to 
be a problem. There are at least three responses to this objection. The first possible 
response is that the interaction between the emerging and energetic (non-relativistic 
and then relativistic) wind with the overlying post-supernova-explosion ejecta may act 
to collimate the outflow. Future multi-d simulations should address this issue in detail. 
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The second potential response is that, in fact, relativistic Poynting-fhrx dominated winds 
can be efficiently collimated, as in the work of Ref. |25j . A final possible response is this: 
given the basic geometry of the magnetocentrifugal wind, it seems natural to suppose 
that the asymptotic radial velocity is largest in the equatorial region. This follows both 
from the fact that centrifugal acceleration is largest at the equator and that the equatorial 
current sheet may facilitate significant dissipation of the magnetic energy. Pulsars (e.g., 
the Crab) provide evidence for high Lorentz factor and energetically-dominant equatorial 
winds |2fil 12 7j . In analogy, is it possible that the geometry of GRBs is "sheet" -like 
(equatorial) rather than jet-like? The solid angle subtended by a sheet with opening 
angle 8 is ~ 6, whereas for a jet it is ~ 2 . This fact has implications for the afterglow 
— a sheet-break rather than a jet-break [2H1 EH — an d for the GRB census. Inferences 
from detailed predictions of the flux evolution may be used to rule out or confirm the 
possibility of a sheet-like geometry in some bursts. I am constructing such models now. 

* * * 
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